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determined using thermogravimetric analyses. Measurements defining 
these events were obtained by comparing scanning calorimeter and 
themogravimetric data at the same beating rate of 10° C per minute. 
The effectiveness and response characteristics of these instruments 
for propellant measurements are also discussed. 

Scanning calorimeter results also indicate that nigroglyeerine 
when mixed in with nitrocellulose does not act independently as an 
ingredients Low temperature measurements over the temperature range 
of the freezing point of nitroglycerine of 1J>° C did not show an 
endothexm 'which ’would be expected if freezing did occur. The fact 
that nitroglycerine begins to evaporate at approximately 165 ° C at 
all pressures confirms this finding. 
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I. PURPOSE 


The purpose of this work is to experimentally investigate reac- 
tions occurring in the surface and subsurface regions of M -2 double - 
base propellants at low pressures (O. J psia to 75 psia). The work is 
performed using the techniques of differential scanning calorimetry 
(DSC) and. thermogravimetric analysis (TSA). It is part of an overall 
program at the University of South Carolina to isolate various critic* 
parameters and determine their effects on the combustion mechanism at 
low pressures. 
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Recently, however, WaescheW and Wenograd(^) have suggested that 
condensed, phase reactions are the controlling factor in combustion 
of composite propellants, particularly at low pressures where Wenograd 
associated .the phenomenon of "flameless combustion" with condensed 
phase reactions- However, the appearance of i visible flams is not 
a requirement for gas phase reactions to occur. Daniels^ 3) 
other evidence in support of condensed phase reactions. He examined 
grains of quenched powder after partial burning and found, by the 
use of dyes, that- there is a changed layer near the surface of the 
powder which is perhaps 0.015 cm deep. 

Rice and (Jinell^ assumed that the powder was appreciably heated 
to this depth and showed how it could not be accounted for by conduc- 
tion from the gas phase. However, they did not define what 
"appreciable" heating was or the effects quenching might have on the 
surface. Evidence for truly condensed phase subsurface reactions, as 
opposed bo heterogeneous surface coupled reactions with the gas phase 
is rather flimsy and as this work indicates, at least for M-2 double- 
base propellants, its significance is very minor. 

Some confusion exists in the literature and among workers in the 

propellants field due to rather loose use of the term "solid phase" 

_ \ 

reactions. This implies to some people that solid ingredients react 
end produce solid products. Simi" arly, the term "condensed phase" 
implies liquid reactants said llqiaid product. Most of the literature 
refers to a decomposition reaction where a solid decomposes to gaseous 
components as a "solid phase" reaction. These gaseous decomposition 


products can then react and produce other gaseous products in the 



so-called "dark zone" and "flame zone." In. the combustion, of double- 
base propellants, it is probable that several different types of 
reactions occur simultaneously. A better understanding of the 
phenomena involved may result if these reactions can be defined 
separately. For the remainder of this thesis, the term "condensed 
phase" reactions will mean no gaseous products axe involved. 

One other term .which should be clearly defined in any discussion 
of solid propellant combustion is "surface." It is unlikely that any 
single plane exists which separates gases from condensed phases. 

Some condensed phase peaks will protrude into the gas side of the 
plane, etc. Also in the vicinity of the "burning surface, " the 
temperature gradients are of the order of thousands of degrees per 
centimeter, and moving a distance of a few microns can change the 
temperature of s plane of "surface" by significant amounts. Thermo- 
couple beads less than 1 mil in diameter are still large with respect 
to these gradients* therefore, locating a plane or "surface" by 
thermocouples is quite difficult. When talking of "surfaces" such 
as the one on which a fflittfesiss surface temperature exists, we must 
remember that this is a plane only and may not be 

relatabla to experimental data. This inability to locate a surface 
naturally confuses the issue of condensed phase heat release if one 


refers to any reaction occurring in the condensed phase as a con- 
densed phase reaction. The amount of heat released in the condensed 
phase will naturally vary with the location of the surface. In a 
high -temperature gradient region where distances of a few microns 


give significant changes in temperature, the reaction rates which 



vary exponentially with temper? ‘/ure can give different answers on the 
amount of heat released up to a given, plane or "surface. " 

Comparison of DSC and TGA data at the same heating rates shows 
that no measurable heat is produced before weight loss occurs. There- 
fore, no condensed phase reactions occur before gas producing 
reactions start- Decomposition reactions producing gaseous products 
are most likely to occur at external surfaces,} however, some may occur 
below the gas- solid, interface at 'surfaces of internal voids or cracks 
in the solid propellant. . The visual observance of bubbles at the 
-burning surface of double-base propellants suggests that some gas- 
producing reactions do occur slightly below the surface. The heat 
measured by the DSC may be from exothermic decomposition reactions 
or from gaseous reactions among the decomposition products which occur 
very close to the surface. Microthermocouple measurements from other 
work at the University of South Carolina, which will be discussed in 


detail later, indicate that reactions do occur yery close to the 
surface. Even though condensed phase reactions do not occur before 
gas-producing exothermic reactions, they still can occur simultan- 
eously with their effects superimposed on the other reactions. 



III'. LITERATURE REVIEW 


Consider the propagation of a one-dimensional combustion wave in 
a premixed gas as being fixed in space with reactant .gases flowing in 
from the left and product gases flowing out, as shown in figure 1. 

The steady- state energy equation describing such a flow as given by 
Frank-Kamenetskii^) for constant thermal properties is 



where 

T = temperature 
x = space coordinate 
A = thermal conductivity . 
m = mass flow rate per unit area 
Cp = heat capacity of gases 
W = volumetric rate of energy release 
If we assume the heat release is descrihahle by an Arrhenius 
expression, then 


W = PQZ expf-EjjRT) (2) 

where R = universal gas constant 

Q = energy of reaction in cal/gm 
P = density 

Z - preexponential factor 
E^ > activation energy 



This premised gas equation has been solved with the following 
boundary conditions: 

T ~ T Q at x = -» T 0 = initial temperature 

T = Too at x = ^ = flame temperature 

The resulting equation, which is given by Fliukhin^ is 



This concept . of flame propagation or deflagration accounts for 
the propagation Of flames through, gases and is considered by some to 
be applicable to any continuum. 

Combustion of double-base propellants is much more complicated 
than the combustion of a premised gas, however, because a solid is 
involved as well as a gas- Even though the double-base propellant : 
is usually considered to be homogeneous, it may not be truly homo- 
geneous, probably being closer to an emulsion of several com- 
pounds because of insolubles in the : propellant and the nature of 
mixing. The decomposition products of' this solid may not be in a 
premixed state- Additional heterogeneous phenomena may also enter. 
Therefore, many other theories involving the coupling of various 
zones or phases of solid propellant combustion have arisen. 

Heller and Gordon, ) using high-speed photography, observed 
that the casbustiod region of burning double-base propellants con- 
sists of three zones. A ’’foam’’ zone less than 1 iniilinieter in depth 
exists right below the surface in «hich the solid transforms to a 
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viscous liquid and “bubbles form. Above the burning surface is a 
large dark zone of supposedly nonreacting gases whose length varies 
with pressure. As pressure increases, the flame gets closer to the 
surface. The third zone, which exists only at pressures greater than 
100 to 200 psia, depending on the particular propellant, is the flame 


In the models of Parr and Crawford' 0 '* and Rice and Ginell'^' 
shown in figure 2, the "fizz zone 7 ' is introduced as an additional one 
between the "foam zone" and "dark zone" of Heller and Gordon. This 

"fizz zone" exists even when there is no flame at low pressures. Both 

' • . • ■■ i 

of these theories applied the conservation of mass and energy equations 
to the experimentally observed structure of the combustion zones. They 
treat each zone consecutively, considering the heat generation in. each 
region and heat feedback from succeeding zones. 

Par? and Crawford neglected the "dark zone" and assumed a first- . 
order "foam zone" reaction and a second-order ’’fizz zone" reaction. 

The fizz reaction was considered to be the controlling reaction at 
lovr pressures where there is no visible flame reaction. They inte- 
grated the combined mass and energy conservation equations for the 
• foam and *\zz zones, matching the temperatures and extents of reaction 
at the interface. A temperature T a was defined to be the temperature 
at the end of the foam zone when no heat feedback from succeeding zones 
occurred. A set of tables was provided for ranges of values of activa- 
tion energy, surface tsrperature of the fizz zone, and extent of 
reaction as functions oir T & /T s . The theory was made to fit burning 


rate versus pressure curves, hoping that the values of activation 



energy, temperature, end extent of reaction thus determined were the 
correct values. The model predicted a limiting very low pressure 
coefficient of 

5 log m 

rt - r- — ° — 

d log P 

As shown /by equation 72 of the cited reference, 

xp = o = O.50 

As the authors themselves point out, the experimental value is O.75. 
Therefore, a good, fit of the theory cannot be expected at low pres- 
sures. It should also be noted that condensed phase reactions (no 
gases involved) are usually considered to be pressure independent, 
that is, zero-order reaction. By assuming first-order reactions in 
the foam zone, Parr and Crawford seem to be implying that the foam 
zone reaction is heterogeneous. 

The approach of Sice and Ginell was quite similar to that of 
Parr and Crawford except that they neglected the "foam zone" and 
used empirical data to define the length of the dark zone at various 
pressures, thereby varying the heat*transf er rate from the llame zone 
as it moves closer to the surface. , The dark zone was considered to 
be a region in which required concentrations for the flame reaction 
/built up- They included a treatment of the fizz zone when no flame 
Occurred. It also is o curve fitting process in which parameters 


m 


(5) ; 


such as surface ten^perature with no feedback from the flame and 



activation energy for the fizz zone were considered to be the values 
required to fit the experimental data of burning rate versus pressure. 

Adams') gives a very interesting treatment of coupling the zones 
involved in double-base combustion. He includes the diffusion equa- 
tions for reacting species with the conservation of mass and energy 
equations* For plane one-dimensional flow* these, can be written as 


m 2a -■ -w. 

dx i 


where M is the total mass flux* is the fractional mass flux 

for component I, x is the space coordinate* and is the net 
rate of loss of component I due to chemical change- The diffusion 
equation is given as 


»i 3 «t “ p % *r 


where p is the mass density and is the mass fraction concen- 
tration of component 1 with maos diffusivity %. The energy 


equation is writ ten as 


- k f =0 


h^j cai^Onent enthalpies per unit mass* are defined by 


\ ~ %i + J 



and is the standard heat of formation of species I. The 

summation in equation (8) is taken over all components of the system, 
and k is defined as the hulk thermal conductivity. 

If we assume no heat loss from the combustion products occurs 
and that the system tends to a state of equilibrium, the derivatives 
.vanish at the hot boundary. The -total: enthalpy' then tends to a 
constant value and equation (8) can be integrated to give the temper- 
ature gradient in terms of the local component mass fluxes and 
enthalpies as 


k li = M X G i b i * M ^ G i h i 


hot boundary 


If the Levis number for each species is unity, then PD^ = k/Cp, and 
if uhe enthalpy tends to a constant value at the hot boundary, 
equations (j) and (8) can be combined and integrated to give the 
result that the local enthalpy remains constant from the cold through 
to the hot boundary. This can be expressed as 


- ^f m 


'-IS 


where £hf m is the standard heat of, formation of the final product 
species, is the temperature at the hot boundary where ■ ^in “■ ^ 

and Cp is assumed to be independent; of composition. This same 
result is derived, by Bice end trinell. ^ 

Adams then applied these equations to a solid A undergoing an 
exotheittic aero-order reaction to give a .gas phase product B which 



is then converted to a stable product D in two stages of first and 
second order, respectively, 


k A 

A — B 
zero 
order 


k- 


■B 


irst 
order 


C - 


seconf’ 
order 




The resulting solution for the mass burning rate is 


-where 


M 8 = 


2k a P a Z a 


Cp(2T^ - Tg - T 0 ) E a 


TT exp(-E A /RT 8 ) 


( 12 ) 


k a = thermal conductivity of solid phase 
P a = density of solid phase 

Z a = frequency factor of solid phase Arrhenius expression 
T s = surface temperature 
T 0 ® initial temperature 

T’ = minimum surface temperature defined by setting the 
temperature gradient in the gas phase adjacent to 
the surface equal to zero (no heat feedback from 
gas phase) 


Adame provides plots of the mass flujtn? for the solid phase. 


and f irst-o rder gas phase reactions as functions of the surface 




The mass flux for the coupled system and the value of T s is deter- 
mined "by the intersection of the two curves. The. solid phase curve 
is independent of pressure, hut the curves for the first-order gas 
reaction are shifted vertically upward as the pressure increases, 
being proportional to the square root of the pressures. Below a 
certain pressure F* the intersection lies on the nearly vertical 
section of the curve where T s -*» Tg. The burning rate is effec- 
tively that of the adiabatic solid phase reaction. Above P* the 
intersections occur on the section of the curve which is nearly 
parallel to the T s axis. The burning rate is then determined by 
the gas phase reaction. He also points out that the exact solution 
of M 2 tends to infinity as T s -> Tg . 

Even though the above theories have been available since 1950 
and experimental measurements to verily their calculated valt 'S of 
activation energy, minimum surface temperature, and heats of reaction 
in the different zones are clearly needed, they have not been avail- 
able. One reason may be a general lack of emphasis on double-base 
propellants in recent years. 

For composite propellants, some work has been performed in this 
- area starting about 1966. Waesche^ and Wenograd*' 2 ^ have attempted 
to determine activation energies for condensed phase reactions. These 
authors assert that the combustion of composite propellants is con- 
trolled by condensed phase reactions. The condensed phase reaction 
is aspumed to follow the Arrhaiius law and the propellant surface 
temperature is aubatituted for the flame temperature, in Arrhenius' 
expression. They assumed aero order kinetics for the solid phase 





and assumed, that the surface temperature was determined by equilibrium! 
vaporization as suggested by Powling and Smith. ' Waesche and 
Wenograd's measurements were taken with a differential scanning 
calorimeter of the same type used in this work. Using the heat 
release data from the scanning calorimeter in the mass burning rate 
equation from the Zeldovitch theory mentioned earlier in this liter- 
ature review, they obtained fair agreement with predicted deflagra- 
tion rates in the low-pressure flameless combustion region. Wenograd 
cited the flameless combustion as evidence for condensed phase 
reactions. However, as mentioned previously, a visible flame is not 
a requirement for gas phase reactions to occur. Microthermocouple 
measurements with double-base propellants at the University of South 
Carolina by iBiompson^®^ ani scanning calorimeter measurements in this 
work indicate that gas phase reactions do occur very close to the 
surface in the case of double-base propellants. Waesche also states 
that' the differential scanning calorimeter is applicable to data 
acquisition for condensed phase reactions "because the placement of 
the sensing elements in the base of the sample holders is such that 
only condensed phase enthalpy changes are measured, and the gaseous 
reaction products are swept away by the steady purge stream. " This 
statement is true only if there sure no gas phase reactions close to 
the surface of the sample iand no heterogeneous attack of the surface 
by the decoaposition products occurs. 

Crawford, Kuggett, and McBradyf^V using a bomb calorimeter 
pressurized with nitrogen and using low loading density of propellant, 
obtained some data on beat of reaction for flameless combustion. 




Their data is reproduced in figure 3- The propellant they used is 
5^ percent nitrocellulose, 43 percent nitroglycerine, and 3 percent 
ethyl centralyte. It provided flameless combustion below 200 psi. 

A heat of reaction of 500 cal/gm was reported for flameless combus- 
tion. They also analyzed the products of r-*;ction and mentioned that 
at low pressures (low in their case meaning 20. '.0 300 psi) an oily 
residue remained in the combustion chamber. This was thought by them 
to be glyoxal. Analysis of the products of reaction also suggested 
to them that nitric oxide and other simple organic molecules react 
in the gas phase very close to the surface at low pressures. They 
felt that approximately one -half of the total heat of reaction is 
liberated in this region and the temperature of the reaction zone 
rises to perhaps 1,500° K. In the final stages of the reaction, 
nitric oxide reacts with the remaining oxidizable material and the 
flame temperature’ reaches its maximum value of approximately 3,000° K 
if the pressure is high enough for the flame to exist. 

In 1366 a translation of a Russian journal article by Zenin^ 12 ^ 
appeared* .His experiments were performed in 1958 using thermocouples 
to determine temperature profiles for deflagration waves. These 
temperature profiles were used to obtain the amount of heat released 
in each of the combustion zones and to estimate the heat supply from 
the gas phase to the condensed phase. The results of Zenin's experi- 
ments are shown in figure 4. His results roughly agree with this 
work even though the methods used in obtaining, them are not given in 
detail and he was hot correcting for thermocouple losses and high 



temperature gradient effects- Such effects can he significant as 
shown by Suh and Tsai.' 1 ^ 

. . i 

Although the heat of reaction and activation energy as a function, 
of pressure are very important pieces of information, especially at 
low pressures due to its rapid change, they are not readily available. 
Lenchitz and Haywood^) measured the heat of reaction of modified 
double -base propellants and found it increased with pressure, approach- 
ing an asymptotic value. Lenchitz also measured the heat of 
reaction of M-2 propellant with and without flame and found the . 
difference between- the two cases to be small. Using a bomb calorimeter 
and low loading density, he obtained values of k39 cal/gm, 588 cal/gm, 
and 783 cal/gm at average burning pressures of 1^-7, 111, and 221 psia* 
Lenchitz 1 3 values agree well with the results presented in this work, 
although the measurement techniques are different. 





the heat contents (B) of the reactants, products, and transition 
state are shown dlagrammatically. The distance AB corresponds to 
the difference in energy "between the transition state and the 
reactants state, which equals the activation energy for the forward 
reaction E^. The activation energy for the "backward reaction Eg 
corresponds to distance BC. The heat of reaction is given "by . 
AC. Then we have 

% - Eg = AH (15) 

Using the Van Hoff isochore (see Liepmann and Boshko^), p. 30) 
which is 


d 2 n K = _£H_ 
' dT rt 2 


(16) 


we subabitute (l 4 ) and (15) in (l6) to get 


Ji 

dT 


{ ^l\ Et - Ec 

r^r~~R~ 


(IT) 


aid rearranging this 


d In % d In kg _ Ej_ Eg 
dT “ “ dT ■ “ rP * Rp 


(18) 


This suggested to Arrhenius that the forward and backward reactions 
inight act independently and he considered the two equations . 


d In k-j Ei d In kn 

-jfe =i —4 + C and . 1 . ..i . j fc 


at 


dT 


a . unm . ff> 0 




where C is a constant. Experiments proved, his suggestion to he 
true) resulting in the simplified form 

d Zn k _ ■ E ' j 

dT rt 2 ^ 

which can he integrated assuming E . is independent of temperature 
to get 


which can he rewritten as 


k a A. exp(-E/RT) 


where A is a constant. This constant is determined by the frequency 
at which molecules with sufficient energy collide and the orientation 
requirements for reaction to take place. The exponential portion 
describes how many molecule s can overcome the energy harrier. 

In deriving this Arrhenius expression, we are using a grossly 
simplified case. In reality, reactions are not so simple sad the 
' stoichiometric equation illustrated hy equation (15) - even if it 
could he written for the complex mixture of reactants and final produ 
of propellant burning - would hot describe the kinetic reactions which 
actually take plane.; The real reactions are more likely to be a com- 
plex sequence of series, parallel, branching, or chain reactions 
involving active intenwediate species which do not even appear in the . 
stoichioj«stri,c equation. ■ 




If we use the Arrhenius expression to describe any part of our 
.reactions, we are describing it as an overall reaction with one 
activation energy. In reality, it is probably a much more complex 
group of reactions taking place as the. experimental results section 
will show. However, as long as our proposed description agrees' with 
experimental fact, it should be considered useful to the extent that 
it does express certain overall aspects of the reaction and allows us 
to predict what might possibly happen under conditions, to some degree 
outside of the range of experimental measurements. Since we are: not 
trying to determine the kinetic details of our propellant reactions, 
but. only the overall effects, the approach' of assuming an Arrhenius 
type describing expression does not seem unjustified. 


t: 
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v. Expmmsmta, apparatus and procedure 


1. Differential Scanning Calorimeter - (DSC) 

A. Description of the. Instrument 
. The scanning calorimeter used In these experiments Is manu- 
factured "by Perkin-Elmer Corporation. The operators manual 
describes; it as follows: "The operation of the. model DSC-1 is based 

on the temperature control of two. miniature sample holders in the 
sample holder assembly as shown in figure 5- ■ A block diagram of the 
instrument Is shown In figure 6. The system consists of two separate 
"control loops, " one for average temperature control, the second for 
differential temperature control. In the average temperature looj , 
a progrsmer provides an electrical signal which is proportional to 
the desired temperature of the sample and reference holders. The 
programer temperature information is also relayed to the recorder 
temperature marker pen and appears as the abscissa scale marking. 

The pfogxamer signal iseompared with the average signal from platinum 
resistance thermometers permanently embedded in the sample and refer- 
ence holders- The resultant difference signal is amplified in the 
average temperature amplifier. If the temperature balled for by the 
progrataer is greater than the average temperature of the sample and 
reference holders, more power will be fed to the heaters of both sampl 
holders -which, like the resistance thermometer a, are embedded in the 
holders.. If the average temperature is higher than that demanded by 
the prqgfjMaef the power to both heaters will he decreased. In this 

ef 13ie holders is made to track the 


"lii the differential temperature control loop, signals repre- 
senting the sample and reference temperatures, measured by the platinui 
thermometers, are fed to a circuit which determines whether the 
reference or sample temperature is greater* The differential tempera- 
ture amplifier, output will then proportion a small increment of power 
between the reference anC sample heaters in such a way as to correct 
any temperature difference between them. This is done by increasing 
the power ' J one while decreasing’ the power to the other. A signal 
proportional to the differential power is also transmitted to the main 
recorder pen. The integral of the resulting peak is the internal 
energy change. The direction of the pen excursion, will depend upon 
whether more power is required in the sample or reference heater. " 

The operator raises the average temperature of the sample holders 
to a desired initial temperature and selects one of eight available 
temperature program rates between 80° C per minute and 0.625° C per 
minute. Most of the run? for these experiments were at 10° G per 
minute or 20° c per minute. When any endothermic or exothermic 
reaction occurs in the sample, the change in power required to main- 
tain the sample holder at the same programed rate as the reference 
holder is recorded. The chart abscissa gives the temperature and the 
area under a curve gives the total energy transfer to or from the 
sample. The anplitude of the deflection from the baseline represents 
the exothermic or endothermic heating rate in miilicalories per 
second. The sensitivity to heating rates of various samples is 
selected from 2, h, 8, 16, or 32 miilicalories per second ranges for 
^ai-7Si8aJ.a ; ' racoirdeS'’- deflection. • . 



Provision -is made for removal , of decomposition gases "by allowing j 
an inert gas to. flow into the sample holder assembly. A sketch of the 
gas flow path is shown in figure 7. Examination of this figure will 
show how gaseous reactions which occur close to- the surface cannot be 
swept sway by the purge gases. The sample holder and its aluminum 
dome-shaped cover constitute an essentially closed system except for 
the small amount of purge gas which leaks in around the crack between 
the holder and cover against the outward flow of decomposition gases. 

A photograph of the; assembled DSC and TGA system is shown in figure 8. 

. A special low-temperature enclosure to surround, the sample holder 
assembly is provided with the scanning calorimeter for measuring heats 
of fusion for materials with low temperature melting points. This 
enclosure is constructed like a wide-mouth vacuum-type ''thermos" 
bottle. Liquid nitrogen, can be poured into this enclosure to surround 
the sides and top of the sample holder assembly with a jacket of 
liquid nitrogen. The sample holder assembly can thus, be cooled to 
-100° c. Programed scanning rates can be initiated at this tempera- 
ture. Crushed dry ice and acetone can also be used as a coolant to 
cool the sample holder assembly to approximately -20° C.' The opera- 
tion of the system is essentially the same with the low-temperature 
enclosure as with the normal enclosure except that extra liquid 
nitrogen must be added periodically.. 

The DSC is calibrated by using a 10 mg lead reference sample. 

This sample is placed in the sample holder and a programed temperature 
rise rate of 20° C per : minute is initiated. -When the melting point of 
lead (600° K) is reached, the recorder begins moving in the endothermic 



direction. If the DSC temperature dial does not read 600° K at the 
start of this endothemn, the temperature calibration control is 
adjusted to reduce the difference and the procedure repeated until 
the dial, reads 600° K. The . area of the resulting endothermic curve 
and the known weight and heat of fusion of lead provide a means of 
calibrating the DSC for accurate measurements of the total heat going 
into or being generated by a sample. This method automatically 
accounts for the resistance to heat flow across the interface between 
the sample holder and the aluminum sample pan which encloses the 
sample. * 

B. Modifications for High Pressure Measurements 

The DSC-1 scanning calorimeter as available from the manu- 
facturer is capable of measurements at atmospheric pressure and below. 
For the work reported in this thesis, pressures above atmospheric 
were desired for comparison with data from other work at the Univer- 
sity of South Carolina* An additional sample holder assembly and 
mounting hardware was obtained from the manufacturer. This extra 
sample holder was movusted in a spherical pressure vessel shown in the 
sketch, of figure 9 and photographs of figure 10. Electrical connec- 
tions are tapped into the existing atmospheric system electronics 
and a switch added to allow use of either sample holder assembly. 
However, to eliminate possible variations, the high pressure vessel 
is used for vacuum and high pressure runs. Purge gas and vacuum 
Connections are connected to. the pressure vessel also: A continuous 

purge flow of 20 to 30 cc/min, as: measured by flowmeters, is main- 
tained during all runs above and below atmospheric. The pressure 


vessel is designed to handle 500 psia safely and was bydrotested to 
750 psi before use. However, at pressures, above 75 psi, it is 
difficult to prevent baseline drift. This is probably due to the 
increase in thermal conductivity of the environmental Inert gas which 
is argon in this case. Argon is used instead of nitrogen because there 
is some ' indication that nitrogen affects the reaction processes of 
doublet-base, propellants. The vacuum, purge gas connections, and valve 
system for the high pressure DSC system are sketched in figure 11. 

C. Sample Preparation aid Experimental Procedure 

Two different types of samples are used in the scanning 
calorimeter experiments. The first type is cut from a rod of pro- 
pellant approximately 1/8 inch in diameter. Disks of 0.025 to 
O.C30 inch thick are cut, using a lathe and knife-edge cutting tool. 
Weights range from 16 to 22 milligrams. Each disk is weighed on an 
analytical balance accurate to 1/10 milligram. Each disk is placed 
in the small aluminum pans provided with the scanning calorimeter and 
covered with a lid of aluminum which is pierced to allow gaseous 
decomposition products to escape. The edges of the pan are then 
crimped with the special tool provided with the calorimeter to com- 
pletely enclose the sample disk of propellant. The sample pan is 
then placed in the sample holder assembly and an empty sample pan 
placed in the reference sample holder.. (This is acceptable in differ- 
ential scanning cialouimetKy where the specific heat of the reference 
material does not have to be similar to the sample as In differential 
thermal analysis. ) The range selection switch is set at either 2 or 
h millicelories par seccaid for full-scale deflection. These are the j 



most sensitive settings of the instrument. She argon flow control 
needle valre is then opened enough to allow for a 20 to 30 cc per , 
minute purge flow rate.. For runs below atmospheric, the vacuum pump 
is started and the bypass valve opened to maintain the desired 
pressure. The argon needle valve is then readjusted for a 20 to 30 cc 
per minute purge flow rate. For runs above atmospheric pressure, the 
valve at the entrance to the pressure vessel is opened and the pressur 
regulator on the argon supply tank is set to the desired pressure. 
Simultaneous adjustment of the argon needle valve and the outlet valve 
in the vacuum line provides a purge flow through the same path as for 
low pressure runs and at the same rate of 20 to JO cc uer minute. 

Af 5er waiting several minutes for the pressure aiid purge flow to 
stabilize, the prcgramer switch is set to 10° G per minute and the 
stare switch turned on to begin the scan. The zero position, control 
is adjusted to position the marking pen in the center of the recorder . 
chart. The slope control is then adjusted to produce a vertical run 
along the center of the paper when no chemical reactions are occurring. 
When exothermic reactions start in the sample, the recorder pen 
deviates frem the centerline and' moves toward the edge of the paper. 
When full-scale position is reached, the run is stopped by cutting off 
the main power switch and argon supply tank valve. Pertinent informa- 
tion is then written, on the chart paper for identification of each run* 
Measurements such as these were made at various pressures between. 

0.3 psia end 75 psia. Results of these runs are discussed later. 


This same technique of sample preparation is also, used for the low 
temperature runs. The sample holder assembly is covered with the ■ 


liquid nitrogen jacket enclosure and liquid nitrogen added. The 
temperature of the samples is monitored until equilibrium is reached 
and then the previous operation procedure is followed. The low temp- 
erature runs are performed at atmospheric pressure. 

The second method of arocple preparation is required to prevent 
self -heating in runs for which the heat of reaction is measured. If 
self-heating of the sanpfU- fly . | id exothermic reactions exceeds the 
amount- required to mail : y - framed heating rate, inaccurate 

results are obtained. Tl: is pro rented by mixing 10 percent or less 

of finely powdered propellant with a material known to be chemically 
inert over the range of temperature and pressure being investigated. 
Polypheny! ether, a liquid used as. a high temperature hydraulic 
fluid and as a high vacuum lubricating oil, is a suitable material. 

It is available from Consolidated Vacuum Corporation^ 1 ®' under the 
tradename Convalex-10. 

. A sample pan of pure polyphenyl ether was scanned over the 
temperature range of 25° C to 300° C at pressures down to 10 mm Hg 
and found to be stable and did not vaporize. A known amount of poly- 
pheny! ether is placed on a glass microscope slide and a pre -weighed 
amount of finely powdered propellant (10 percent or less ) is mixed 
with the 'polyphenyl ether and stirred to form a slurry. Wrap several 
turns of Mo. 3Q gage nickel wire around a 20d common nail and cut off 
one turn of wire to form a ring which will just fit into the aluminum 
sample pan. A email droplet of slurry is then placed inside of the 
ring, covered with a lid pierced 10 times with a hypodermic needle 
and crimped at theedges. The wire ring prevents the liquid from beii 



squeezed out during the crimping process. The sample pan, ring, and 
lid are weighed "before adding the slurry and afterwards to determine 
the amount present in each sample pan accurate to l/lO milligram. 

Slurry weights of roughly 20 milligrams at 5 to 10 percent propellant 
produce .curves which remain on scale for range settings of 8 to 
16 milliealories per second, depending on the pressure. The procedure 
for setting of pressure and purge flow is the same as for dish samples.! 

• I 

For these runs the programed temperature rise rate is allowed to con- 
tinue until exothermic reactions cease and the marker pen returns to 
the baseline. The "baseline sometimes drifts slightly in the exothermic 
or endothermic direction, depending on how accurately the slope con- 
trol is adjusted at the beginning of a run. Several nans at one 
pressure usually allow a proper setting to be achieved. Some readjust 
ment is usually required when the pressure is changed between runs. 

The heat of reaction runs usually produce a bell-shaped curve. 

■ 

The area, under the curve bounded by a baseline drawn tangent .to. the 
lower edges of the bell is measured with a mechanical planimeter. 

The planimeter is then used to measure the area of a rectangle which 
is one-half of the. chart in width, and in length equaling 1 minute 
of run time. (For example, if the chart speed is 4 inches per minute, 
the rectangle is 4 inches long.) A direct measure of the beat of 
reaction • in calories is then given by 

Calories - (Area of seagle peak)(renge settlng)(6o) ( 2£ y 

(Area of rectangle) (1000) 

This is then divided by the weight of propellant in the sample to get 
the heat of reaction in calories per gram. 
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The procedure given in reference 19 was also tried. The propellant 
■was dissolved in acetone and coated onto an inert solid* cab-o-sil 
(obtainable from the Cabot Corporation* Boston* Massachusetts). For a 
known initial amount of cab-o-sil on which a 10 percent solution of 
propellant. ,in acetone was poured, the weight increase was never 10 . per- 
cent. Similar mixes produced variable weight changes and nonrepro- 
ducible heats of reaction. Orders of magnitude of measurements were 
similar to the polyphenyl ether results, but no definite trends were 
recognizable. 

• One other procedure* using small amounts of pure powdered M-2 pro- 
pellant (less than £ mg), was tried- The resulting curves will remain 
on scale for range settings of 32 millicalories per second. These 
samples also must be enclosed in sample pans with pierced lids to pre- 
vent powder fragments from getting out of the sample pan due to purge 
gas flow or decomposition gas flow and causing spurious peaks on the 
recorder. 

2. Thennogravimetri c Analysis System - (TGA) 

A. Description of the Instrument 

A Ferkin -Elmer Model TGS-1 thermogravimetrie analysis system 
’ is used for weight loss versus temperature measurements. It provides 
a . record of microgram level weight changes in a sample as a function 

a 

at temperature from ambient to 1*000° C. It utilizes s, Cahn "BG n 
electrobalance which is mounted in & glass vacuum chamber* permitting 
. control of the atmosphere around the sample suspended from the bal- 
ane® > ..m. of the elect rebalance into a furnace below it. The Cahn 


eieetrobalaace provides a sensitivity of 0.Q001 mg. The furnace 



outside of a ceramic cylinder 3/8 inch in diameter. The sample is 
suspended in a sample pan inside the heated ceramic cylinder. The 
location of the furnace within the chamber permits close coupling of 
heat source and sample to permit rapid scanning without information 

loss. The low thermal mass of the furnace, allows rapid programed 

' * / o 

temperature rise rates at 11 selected linear rates from 0.3 to 320 C 

per minute. The furnace acts as a heater- and as a temperature sensor. 
In- the temperature sensing mode , ' it forms one side of a bridge circuit 
The other side of the bridge circuit is driven by the outp ; signal 
from the temperature programer. An error signal proportional to the 
tenperature error is developed and fed into an amplifier. In the .. 
heating mode, the amplifier output is connected to the heater by 
means of an electronic system which provides 60-cycle power pulses | 
designed to correct the temperature error. j 

Temperature calibration is accomplished by heating ferro- 
magnetic standards, held in a magnetic field, through their Curie 
temperatures. The magnetic field is provided by a permanent magnet 
mounted to a swinging beam which is moved close to the furnace 
assembly during calibration. The resulting change in magnetic force 
as the Curie point temperature is scanned appears as a weight loss 
on the recorder at the corresponding point in the temperature scan. 
Controls of the TGS -1 are used to adjust the power to the furnace 
so that the indicated, teaperatures of the programer approximate the 
hnown Curie point t«pera , bHires. 

The Cahn electrobalBnce is calibrated by placing Impwn weights 
in pans suspended from the beam and adjusting l^e system controls 



with and without a sample pan until the balance is properly zeroed. 

The weight of a sample is determined from the recorder trace and 
.control settings, depending on Idle selected range of weight change 
required for full-scale deflection of the recorder. Nominal sample 
size is 0.5 to 10 milligrams. Mass range settings of 1 milligram to 
200 milligrams with a recorder range of 20 micrograms to 20 milli- 
grams full scale are available. Argon is used as the inert environ- 
mental gas for the measurements in this work. A photograph of the 
system is shown in figure 8 and the argon gas system connections sire 
in figure 11. 

B. Sasrole Preperation arid Experimental Procedure 

Procedures are based on the operators manual provided 
by the manufacturer. A prime factor in sample preparation is to main- 
tain tbi best possible thermal contact with the sample pan to give 
sharp and reproducible results. Sample sizes are used which fall In 
the optimized range of the instrument of 0. 5 to 10 milligrams. The 
runs at 10° C per minute used a sample of 7 to 8 milligrams. These 
were cut from the same rods of propellant used for the BSC sample 
disks. These were cut to a thickness of about 0.010 inch to provide 
as nriich surface contact with the sample pan as possible. Aluminum 

sample pirns, approximate Jy the sea e size ss BSC sample pans, were 

. ■ . * *• ' » 

first used but could be used cniy once. The carbonaceous residue, 
which ramained after & rim in which complete deflagration occurs, is 
difficult to remove without: banding or damaging the soft thin pans. 

The platinum sss^uLe panswere found to be more sy, » fbr re-use. 



since they could Toe easily cleaned for re-use by scraping with a 
small screwdriver point and washing with acetone. 

For the low scan rate runs# the thin disk samples of approxi- 
mately 7 or. 8 milligrams are cut down in size by cutting off the 
edges to make a snail square sample of slightly over 1 milligram 
weight. Using a full-scale deflection setting of 1 milligram with 
this size sample produces a sigmoid- shaped curve. Several size 
samples and ranges were tried in order to get the residue weight to 
remain on scale sad utilize the full chart width for good resolution. 

The operating procedure is, as follows: We assume that the proper 

tare weights, required for the type pan being used (aluminum or 
platinum), have bean placed on ; the' balance - and the preliminary gross 
balance calibrations made. With the mass dial at 0. 0000 and the 
recorder range dial at 1, zero the recorder using -the set 0/10 
control. Set . the mass dial at 0*5000 snd place a 5"milligram pre- 
cision weight in Hie sample pan and zero the recorder with the set 
5 control. Repeat ’this process with and without the 5-milligram 
weight until no change in the zero position occurs (without adjustment 
of either the 0/10 or set 5 controls). Then with the weight in place, 
the recorder range set to 1, and the recorder still at zero, rotate 
the mass dial one full turn clockwise to move the recorder pan to 
approximate full-scale deflection position. Bet the recorder to read 
exacl^y fUJil scale with the calibrate recorder control. After this 
calibration, the value of the weight on the sample pan will be given 
by the mass dial setting ainusthe recorder feeding to the right of 
position- ' 



j A sample of . the desired size, as determined lay pre-weighing on . 

ah analytical balance accurate to l/lO milligram, is placed in the 
sample pan. No lid is used on the sample pan as in the DSC. A lid 
is required in the TGA only when sputtering of the sample may occur, 
causing some of the sample to fall from the pan. The range is now 
set so that a 1-milligram weight loss will cause a fall- scale 
recorder deflection. This same range is used for all sample sizes in 
these, experiments and for all scan .speeds used. A recorder chert 
speed proportional to the temperature scan speed which will provide 
a trace without extremely small or extremely large slopes for ease 
of data interpretation is selected. The scan rate is selected and the 
. switch turned on to start the heating program. For runs at low heat- 
ing rates of 0.62 titd 1,85° C per minute, a limit control to stop -the 
programed rise rate and automatically cool it bach to ambient is 
available. It is set for a temperature high enough to cover the 
. desired range. This is convenient for long, time -consuming runs, 
which can he set up end left overnight to be completed automatically. 
The resulting data is treated in the results and sample calculations 
sections. • ' 



VI. ESEERIMBHTAL RESULTS 


1. Differential Seaming Calorimeter Results 

, A. Vaporisation of latro^i^erine from Propellant 

In the maxSawsa rune using the 20-iailligram disk 

samples, the first recognisable departure from the baseline occurs at 


approximately lw° C for a beating rate of 10" C per minute. ‘Shis exd* 
thermic indication continues to increase gradually as shown in figure IS 
At approximately 165 C, the established trend in. the exothermic direct 
tied then dips slightly and then resumes its increase . This appears tc 
be. an endothsrm on top of the established exothermic rate. If the mal 
power switch is cut off prior to l65°Q t the propellant sample holder | 
Is found to be clean. But, if the temperature is allowed to increase 
beyond this point, a small amount of brown oily residue is found on t 
imder side of the aluminum dome-shaped cover over the sample holder cu 


and around the outside of the 


holder cup . This phenomenon con- 


oistentiy occurs at 1 65° C for ail. pressures from 0.5 paia to 75 paia 
at a heating rate of 10° 0 per miniate. At this same temperature and 
beating rate , TOil remits (which -will be discussed in detail later) 
indicate a sudden: increase in the rate ofweight los% supporting the 
■ idea of raporitation ■ starting , 'at thi s point . ■ Vov ’.funs at 5? 0 per 
■minute- the anomaly occurs at 157° <3> and at 20° 0 per, minute it occurs 
.at,i73f 'This vaflatio^ wlth heating rata ls a ehareurteristic of : 
reactions .with hi^ aCtivatlea energies (clow reaction rates) which - 


will be 


'in' a later- aeeilea* 


®a> perfojmren infrrji'M. s^aayaiB, & snsple of this brown oily ' 



the sample holder assembly. An infrared spectrum of this, liquid 
indicates it is nitroglycerine, as shown by figures 13 and 14. The 
appearance of this brown oily- liquid at low pressures has been noticed 
by Clary and by Crawford, Hugget, and McBrady. ^ 1:L ^ It is thought 
to be glyoxal by Crawford, Hugget, and McBrady- The infrared spectrum 
of glyoxal is shown in figure 15. Glyoxal is fluorescent under ultra- 
violet light. A sample of the brown oily liquid, if exposed to ultra- 
violet light, does not fluoresce. 

After noticing that M-2 propellant with approximately 2Q percent 
. nitroglycerine produces this perturbation at I65 0 C, it would be 
expected that JEN propellant with 40 percent nitroglycerine would 
produce a greater effect, and propellant No. 691 with 60 percent 
nitroglycerine even more. The JEN did not follow a similar path as 
shown in figure 16, hut the No. 691 propellant shows signs of 
different behavior beginning at about 162° C» It was noticed, how- 
ever, that the, amount of brown oily residue- remaining in the sample, 
holder assembly increased with JIN and even more for No. 691. The 
behavior of a pure nitrocellulose s&nple is also shewn in figure 16^ 

No brown oily residue at all was found in the ssssple holder assembly 
after thepure nitrocellulose test. 

Since the onset of significant heating begins after this apparent 
vaporisation, it is possible that hew reactions are starting- The 
M-2, JE5N, and No. &91 propellaate sll have different ingredients and '' 
different pareentages of MtrocellUlose sad nit3«glycefine. the 
rffftCtic®a oe<lurring iaffly therefore be dlffersit. Scmve may b© 



endothermic and some exothermic. They may also be affected by the 
physical arrangement of these different Ingredients in the nitro- 
cellulose™nltroglycerine matrix. A distinct endotherm which would 
be expected from pure nitroglycerine would not be likely to occur 
wb'in the molecules are distributed id the propellant and mechanically 
constrained. The perturbation we are talking about is also quite 
small and only detectable when using a larger than normal sample at 
the most sensitive heat detection' setting* The behavior of the JPN 
is probably a more normal behavior. The detection of the very small 
perturbations of the M-2 and Mo. 69 I propellant a was indeed fortunate. 

B. Heat of Reaction Versus Pressure 

Typical recorder outputs for the heat of reaction runs are 
shown in figure 17. The heat of reaction varies with pressure as 
. shown in figure 18. Jumps in heat of: reaction occur between Ik . 7 psia 
and 30 psia, between 1 psia and 2 psia, and again between 110 psia 
and 220 psia from Lenchifcz's data points. The results using pure 
propellant powder are shown by X's on figure 18. Possible self- 
heating effects and the difficulty of controlling heat input and 
distribution in the powder, as opposed to the powder in polyphenyl 
ether, probabJy makes the powder results less precise. They are 
included to show essential agreement with the polypheny! ether results. 
This indicates euoy reaction between the propellant and/or its products 
1 of decoasposttion wlth the polypheny! ether,!? it exists at alV is 
very minor.' 3M a Jump j^eaoasena agrees with the results of Suh and 
Glaseyf^J, shown in fi45dre 19> aad the results of Thoeqpsm and Suh» 

The slopes oftlisa curves in figure 19 are parallel froa S.k to 




14-7 psia with a slope change between 1^.7 aid 65 psla. Prom 65 psia 
to 115 psia, the slopes are again roughly parallel and another change 
in slope appears between 115 psia and 165 psia* The data of Thompson 
and Suh indicate that two temperature jumps occur very close to the 
surface and move toward the surface as pressure increases, as indicat 
in figure 20. The pressure at which the jumps in temperature touch 
the surface agree with the jumps in heat of reaction and the change in 
slope of the burning rate versus initial temperature curves. The plot 
of burning rate versus pressure on figure 18 indicates the possibility 
of a jump in the heat of reaction and burning rate at lower pressure. 
However, this one has not been -verified by thermocouple data and is 
therefore speculative. 

The curves of figure 17 are also used to calculate the activation 
energies as a function of pressure given in table I. The data from 
the scanning calorimeter is in the form of distances between the heat 
of reaction curve and a baseline : at the associated absolute tempera- 
ture. This distance is proportional to the rate of heat release from 
the sample and is therefore proportional to the rate constant. Pro- 
vided that a plot of the log of distance versus l/T is a straight 
line, the activation energy can be rpproKimately calculated freon the 
expression derived in the sample calculations section. 


■ indi-Indg . ( 

■■■■i "T ■ '■■■■■' 

where . end dg aredi stances fre® the baseline to the; heat, of ' 
reaction curve at t«af^ratures . %. . asad ; Tg. .iecsraaey,; /' 



values at the extreme ends of the Straight line portion are selected 
to substitute in the above equation. A straight line usually results 
from points at first detectable deviation up to the point where the 
rate of rise starts its inflection to go over the top of the bell- 
shaped curve. 


The derivation of this equation does not make allowance for the 
mass changing. Therefore, in addition to the approximate nature of 



of error. A more detailed discussion of the changing mass effect is 
given in the section on the computer program for DSC simulation of 
the appendix. This is a basic question, which up to now has not been 
pointed out in the aveilablt. literature on the c. lication of the DSC 
to solid propellant decomposition studies. Forttme'-Fly, the results 
of the computer simulation for the DSC heat of reaction curves, 
including the mass effects, indicate that the results of using this 
equation . are fairly accurate. 

C. Nitroglycerine Freezing Point Measurements 

A 20-milligram disk prepared as described previously is placed 
in the freezer compartment of a household refrigerator overnight. 

This sample is removed and immediately placed in the sample bolder 
assembly. Hie low temperature enclosure is then placed around the 
sample holder assembly. A mixtere' of dry ice and acetone Is put in 
the low temperature enclosure. After holding the temperature of the 
sample holder assembly at -20° C for approximately 1 hour, a 20° C 
per minate scan is initiated. A 20° C per minute scan rate allows a 



larger temperate ire range to "be covered without adding more dry ice to 
the low-temperature enclosure and upsetting equilibrium conditions. 
When the freezing point of nitroglycerine at 13° C is passed, no endo- 
thermic. changes are seen, even with the calorimeter, set for maximum 
sensitivity. With 20 percent nitroglycerine in a 20-milligram sample, 
any freezing should have been easily detected. This procedure, when 
repeated using liquid nitrogen as the coolant -and holding . the sample 
at -100° C for 1 hour, produces "the same result (i.e. , no indication 
of nitroglycerine freezing ) . 

It is possible that the nitroglycerine molecules influence the. 
nitrocellulose matrix in a manner similar to the way in which lubri- 
cants affect some organic plastics. The plastics become soft just as 
propellants become softer and more pliable as the nitroglycerine con-' 
tent is increased. JPN with nearly 40 percent nitroglycerine is much 
softer than. M-2 propellant with 20 percent. Regardless of the 
mechanism involved in the softening of nitrocellulose by nitroglycer- 
ine, it appears that the nitroglycerine does not exhibit the same 
properties as a propellant Ingredient- as it does when it is in pure 
form. 

The liquid nitrogen cooled experiment, if allowed to continue, 
shows first signs of exothermic reaction beginning between 60° C and 
70° C, as shown in figure 21. This agrees with the first signs of 
weight loss on TGA runs at the same heating, rate, and indicates 
that only surface reactions are involved. If solid phase reactions 
without gasification existed, heat would be released before weight 
loss, but this does not occur in any easel The temperature at which 



weight loss begins (for the same heating rate) is at least equal to 
or slightly less than the temperature at which exothermic heat release 
is detectable. The low temperature runs are used to detect the first 
signs of heat release. If the normal enclosure for the samples 1 b 
used, sufficient time is not always available to finish adjustment of 
the slope control at the maximum sensitivity settings and 10° C per 
minute rise rate. Starting out at the lower temperature allows plenty 
of time to establish a straight baseline before reactions occur. 

2. Thermogravimetric Analysis System Results 

A. Vaporization Of Nitroglycerine from Propellant 

With an approximately 8 -milligram propellant sample at a scan 
speed of 10° C per minute, the first detectable weight loss begins at 
approximately 60° C. The weight loss rate ' increases slightly with 
temper 8,ture as shown in figure 22 up to 165 ° C where a sudden increase 
in rate of weight loss Occurs. This increase in weight loso rate 
coincides with the endothermic anomaly from the scanning calorimeter 
results. This change in slope occurs at the same temperature for 

different pressures. This is another indication of a possible vapor- 

' . . • ' . r 

ization beginning* or at least, a new type of reaction associated with 
the vaporization of nitroglycerine. The fact that the temperature at 
which this occurs does not change with pressure is again an Indication 
that the nitroglycerine does not act independently as an ingredient 
wien it is mixed in the propellant.. However, since the particular TGA 
systemcan only be used at atmospheric pressure and below, no pres- 
sures above atmospheric were investigated. t 
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After the "knee" is reached, the weight loss rate increases very 
rapidly, and rapid deflagration of the sample results. Cutting off 
the main power switch after the "knee" does not prevent complete 
deflagration from occurring as could be done with the seaming 
calorimeter. Evidently the larger thermal mass of the TGA heating 
system keeps the temperature up long enough for runaway conditions 
to occur, whereas the low thermal mass of the scanning calorimeter 
cools rapidly enough to prevent it. 

B. Determination of Order of Be act ion 

If a small sample of propellant (l milligram J is placed in 
the TGA sample pan and exposed to a very slow heating rate of 1.25° C 
per minute or less, the sample does not reach, the runaway condition 
of deflagration. Instead, a curve shown in figure 2j results. A 
solid carbonaceous residue of approximately 10 to 15 percent of the 
original weight remains. It is believed that this carbonaceous 
material could also be a product , in the flame zone of a burning pro- 
pellant , since glowing carbon particles axe what cause the flame to 
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be visible. Ucing the method of Kusng-Hua , a plot of 


A to — 
dt 

A to SL 


to — 

versus -—3L. 


produces «. straight line and the order of reaction is given by the 

A ' dJf 

' A' to Jg " 

intercept on the - ™ -- axis. 

A to W r 

The derivation of this eguaticu is r^rdduced in the sample calcu- 


lations secticH' . the order of reaction appears to follow an 


if^sn 



apparent zero-order .trend lip to. 16;'" Z where it . then . changes to 
approximately second order as shown by figure 2k* This curve is 
based on slide rule calculations.. T:. . rdsr to achieve more accuracy, 
since some of the mass ratios involved are very close to 1 in value 
and close together, a computer curve fitting method is used. A 
description of this program and the Fortran statements are included 
in the sample calculations section. 

• The computer results indicate a zero -order reaction up to the 

o 

"knee" of- the curve, at approximately 150 C and then a change to 
approximately a second-order reaction. Portions of the computer 
printout sheets for this, curve are shown in figures 25 and 26. 

Figure 25 is in the vicinity of the "knee" of the cur^e, which occurs 
at a weight value of approximately 0.98k milligram in figure 25* 
Values in column l are the assumed orders of reaction which the com- 
puter is told to ex amin e for heat fit to the experimental curve. 
Columns 2 and 5 are the weights at the ends of the interval under 
consideration. Column k is the activation energy which makes the 
assumed order fit the experimental data best. Column 5 is the pre*> 
exponential factor for the best fit. Columns 6 and 7 are the actual 
error and a root mean square absolute value error. Column 6 is the 
one the eomputer uses to select the best overall fit wM.ch is given 
in the last row for a given interval (after the highest assumed order 
row). An overall best fit is not always calculated because of the 
characteristics of the programs unless column. 6 shows a positive and 
negative value within the range of assumed Orders. Hotlce that no 
best fit value is found for the interval 1.03k to 0-984 of for 



intervals prior to that, but the least error is at zero order and 
decreasing- Over the next .interval, the order goes to 1.687 and 
increases thereafter. The average value of the order for the next 
10 intervals is 2.162 and the sw rage activation energy over this 
interval is , $00 cal/mole . At the end of these 10 intervals, the 
weight of the sample is down to 0 . 454 - milligram and the temperature 
is approximately l8o° C. The intervals ■ are equally spaced, based on 
0.050 milligram weight loss per interval,. The change of order at 
the "knee" of the curve is not unexpected, based on the premise of 
nitroglycerine vaporization beginning at this point as with the other 
TGA and DSO curves. The occurrence of this "knee" near I50 0 C as 
opposed to the "knee" in figure 20 at 165° C is due to the lower 
heating rate of 1*25° C per minute. 

The same computer procedure was attempted on the samples which 
reached runaway deflagration conditions at the higher heating rates. 
Hue to the rapid changes of weight loss after the "knee" of such 
r curves, good inputs for the computer, which works over l ill intervals; 
| are hard to obtain from the curves. However, based on the trend of 

error decrease, the order was nearer zero before the "knee" and 
changed to higher order after it. 

5* Condensed Phase Heat of Reaction 

The original intent of this investigation was to determine- the 
effects ofpressureoncondensed phase reactions, The use of the 
differential scanning ealoriaster, which has been used by others to 
obtain such ^ information, is apparently inadequate, at le&yt for double” 
base propellents.. The 3n®^ido then fidrises as to hotr one can obtain. 



the condensed phase heat of reaction or heat of decomposition. This 
can he obtained by the following Interpretation of autoignition meas- 
. urements made at the University of South Carolina as described by 
Suh, Tsai, Thompson., and Moore. 

Cylindrical rod3' of'M-2 propellant like the ones used, to make the 
disk samples described in the experimental procedures section of this 
paper were coated with Bpon Adhesive 9*1-6 (0. 005 cm thick) for planar 
burning. The coated propellant was cut into 2.75-inch lengths. Into 
each length two No. 80 size holes for indium burn wires and one No. 75 
size hole for the thermocouple were drilled. The propellant was then 
placed in a stainless steel test chamber and connections made, as 
shown. in figure 27. The indium wire used was Ind&lloy Intermediate 
Solder No, 6, 50 gage, with melting point Of 280 to 285° C. The 
thermocouple was made of 30 gage, single stranded, enameled iron and 
constahtan wires joined together by resistance welding. 

After having completed the above steps, the propellant temperature 
was raised slowly to some desired value by passing heated air around 
it, and the temperature history was recorded. Once the propellant 
temperature reached the desired level, it was ignited by a hot wire 
heater made of platinum wire. As, the flame front passed through the 
sections containing the indium wires, they melted, yielding the 
information needed to calculate the burning rate. The above procedure 
was. used with air as the surrounding atmosphere for the propellants 
For the series of tests done with argon as the surrounding atmosphere^ 
2 minutes before igniting the propellsnt thc air supply was turned off 
tsnd the argon supply was turned on.. . 



Figure 28 shows the experimental results for burning rate as a 
function of initial temperature. The experimental results obtained 
using both azvpn and air fall on the same line. The asjraptotic 
temperature at which the burning rate approaches infinity is found to 
be 1^5° C ( 290 ° F ) . This is defined to be, the autbignitioh tenperatur 
This experiment provides us with a very useful piece of informa- 
tion, which is the minimum surface tenperature at the solid surface 
with no heat feedback from the gas phase. The experimental arrangement 
indicates that if a piece of propellant is heated to 1^5° c, it is 
just about to ignite by Itself and the slightest amount of energy 
which is added to it externally results, in sudden ignition and very 
rapid deflagration at a rate approaching infinity. Theoretically, a 
burning surface temperature of 146 ° C would be sufficient to insure 
burning of a propellant preheated to 145° C at a rate approaching 
infinity. This experimentally determined value can be used in the 
theory of Adams^) to obtain the pure condensed phase heat of reaction 
Hex-e we rewrite, for convenience, Adams' equation from the literature 
review section, that is, 




KprBMdJjr, is the stetndaspd heat of formation of the final product 

species, and is the flapte i^^erdture at the hot Irouhdhry. If 
represent s the heat of formation of reactant ^eeies, the left 
side of efoation^ As the definificm of heat of reaction. Since 
at the. hot bhuEdary « 1 we can tafee m average value of Cg, 



considered as a constant,, and Cp(T m - T) gives us the overall heat 
of reaction for the propellent burning. However, if we are interested 
in the condensed phase heat of reaction, only, we now have a value for 
which the temperature gradient at the solid surface is zero, com- 
sponding to the temperature gradient of zero at the hot. boundary 
in the above expression. . This value is the minimum surface tempera- 
ture value obtained by the autoignition measurements. Equation (ll) 
can be applied to the condensed phase region considering as 

the heat of formaticai of condensed phase products and . as the 

reactant species again. How our hot boundary temperature is lt-5° C 
and Cp(i^5 - .T ) gxves us the heat of reaction for the condensed 
phase reaction (i.e., no gaseous products formed). , 

If we wish to determine the heat of decomposition for the reaction 
which produces gaseous decomposition products, we would have to deter- 
mine the temperature of the gases which result from the decomposition 
before these products react to form the final gaseous product species. 
If this temperature were known (without any heat feedbacfc from gaseous 
product reactions ), we could follow a similar procedure as before to 
_ deteiiSirjte the heat of decesjpb sit ion for the reaction producing gases, 
from the solid propellents As yet this value is not available. . 

" Qiiie'v f&r M-2 propellant sds an 

: it^tiai of 2*j° C is therefore given' by 

H a * b. 5?(l^5 ~. 25) « 45 eal/gpij . (2V) 



It is thus a single matter to determine the condensed phase heat of 
reaction for any propellant for ''which the autoignitioh temperature 
and specific heat can he obtained. 

Adams' theoretical equation for the mass turning rate discussed 
in the literature review section, that is. 


M 2 = ^aPaZa 


Cp( 2 T s - Tl - T 0 > E a 


■fft) 


( 12 ) 


also indicates that if M-2 propellant is preheated to the minimum 
surface teaperature, then I s = T 0 = Tg . - 145° C, and the denominator 
of equation (.12.) goes to zero end an infinite turning rate results 
Equation (5-6 ) of PliuJihin^ which is 


.. ■■ g^jws . Xfej . 

E s(S T s ~ V'o ) 2 CGpT s - Cpf 0 ) 2 


( 25 ) 


also inpiies this tj^pe tehav lor. Since dT a /dx equals zero at the 
surface under no heat feedback conditions, the second term equals 
zero. if 2 0 .» in the denominator of the first term, the mass 


t?rmiag rate again goes to infinity. 

®iis value of appro^twately 45 cal/gm agrees in order of magnitude 
with Senin's^ 2 ' value df QO cel/^. The composition of his pro- 
pellant is not given, and his results are based on thermocouple data 
uacorrected, for leadlosaes. '2enih'' s' method of ohtainiiig the heat 
release below the surface from thermocouple d&ta would include heat 


from pure condensed 


reactions and from any hetero- 


geneous orgasphaae reaction a occurring aivolds orother surfaces-;' 




"below the solid gas interface. This would naturally "be expected to 
give a higher value. Zanin also indicates that if the pressure is 

■ • t 

lowered enough ( 0 . 5 mm Hg) gas effects are negligible. If the 
line of pressure versus heat of reaction in figure 18, is extrapolated 
to 0.5 nnn Kg, a value between 50 and 100 calories per gram is obtained. 
This is given as a matter of interest only, because extrapolation over 
three orders of magnitude is rather dubious. 




VII. COMPUTER SIMULATION OF DSC 


As discussed previously in the DSC results section, there is a 
time dependence associated with the sample under test for different 
heating rates. The apparent endothermic auGtnaly relating to the 
vaporization of nitroglycerine which occurs at 1 65° C with a heating 
rate of 10° C per minute shifts to 157° C for a heating rate of 5° C 
per minute and to 173 P C for a heating rate of 20° C per minute. The 
hell- shaped heat of reaction curves also shifted at different heating 
' rates, as shown in figure 30. 

In order’ to determine if this is due to the time dependence of 
the sample rather thear a characteristic of the DSC, an equation was 
written describing the heat release rate sis a function of time 
including mass change effects. This equation, for which the deriv- 
ation and method of solution are given in the sample calculations 
section of the appendix, is 


Q=/ — .. exp (- exp (- -^tdT 

Vo i + T 2 ^ v.«y \ **/ 


( 26 ) 


, The digital, computer results for the solution of this equation are 
given in figure 31. The changes in the locations of the peaks and the 
shapes of these curves approximately agree with the experimentally 
detesBiined ckeews of figure 30^ 

The DSC value for activation energy of 51,900 cai/aole and fre- . 


queacjr. factor of lO^ did not produce the desired curves. It was 
neeessfiry to reduce the activation energy to 45,000 cal/mole and the 

in order to- get the curve to, -return close to 




zero on the high temperature side of the hell shape. The corresponding 
values for the mass loss expression which produced the desired curve 
shapes were = 56 , 000 cal/mole and A = 10^ which roughly agree 
with the experimental values as shown in figure 26. 

The fact that the independently obtained DSC and TGA experimental 
values agree this well and could he used together to reproduce the 
heat of reaction curve shapes provides some assurance of their 
accuracy. The use of the Arrhenius expression and equation (A~5) 
seems justified , and the shift of the peaks is due to the propellant 
and not the DSC. 

This procedure provides a method of checking activation energies 
which accounts for mass loss effects and is probably more accurate 
than present methods. 
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VIII. CONCLUSIONS 


Experimental measurements of the overall heat of reaction for 
M-2 double-base propellants at low pressures are presented for the 
first time. Agreement with independently obtained results of burning 
rate and thermocouple data of others indicates the values are 
reasonably accurate. 

The differential scanning calorimeter, when used independently, 
cannot measure condensed phase heat generation in double-base pro- 
pellants because of gas phase reactions which occur close to the 
burning surface. It can be used in conjunction with thermogravimetric 
analysis measurements (taken at the same heating rate) and high 
response microthermocouple measurements near the burning surface. 
Comparison of data from these three sources indicates that the main 
heat producing reactions in M-2 double-base propellants are gas phase 
reactions occurring near the surface and heterogeneous surface 
reactions which may occur at external surfaces or at surfaces of 
voids and cracks below the burning surface. 

Vaporization of nitroglycerine from the propellant has been 
verified. The beginning of vaporization coincides with the beginning 
of significant heat release and weight loss rate. The probable order 
of reaction is zero up to the point where vaporization begins and 
then a change to approximately second order occurs. 

An indirect method of determining the condensed phase heat of 
reaction is given. It is based on the determination of an arti- 


ficially defined minimum surface temperature using autoignition 


measurement s . 



IX. RECOMMEHDATIONS 


The next area of investigation should directed toward measuring 
the amount of gas-producing reactions occurring below the gas- solid 
interface and possibly separating this from the condensed phase con- 
tribution. This has already been initiated at the University of South 
Carolina. 

A very low ‘pressure chamber for use with the differential scanning 
calorimeter could be constructed to measure the heat produced at 
pressures of 0-5 mm Hg and below. If gas phase reactions can be 
eliminated at such pressures as claimed by Aleksandrov ^? 1 ) # the heat 
of reaction may reach a constant value. However, it could be just 
the reactions among gaseous decomposition products which would cease 
at these pressures. Decomposition reactions may still occur and 
produce gaseous products which do not react further. If a TGA chamber 
for these low pressures could also be constructed, weight loss meas- 
urements in conjunction with the DSC measurements at very low pressures 
may determine if this is occurring. The true condensed phase heat of 
reaction can only be obtained if heat changes can be measured while 
no gases are being produced, that is, no weight loss occurring as 
determined by TGA measurements at the same pressure and heating rate. 

High pressure TGA measurements should also be performed. The 
jumps in burning rate shown on figure 17 should be indicated by a 
change in weight loss rate in this pressure range. 

Individual ingredients used in the composition of double-base 
propellants do not act in the mix as they do independently. Therefore, 



in future studies less emphasis should he placed on individual 
ingredient studies and more on the mixfaire. 

An overall program similar to the one which has been performed 
on double-base propellants at the University of South Carolina should 
be repeated for composite propellants. This should include combined 
DSC, ISA, and high response microthermocouple measurements as well as 
autoignition and burning rate versus initial temperature measurements. 
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Figure £>.- Models of Parr and Crawford (a) and Hice and Ginell (t>) . 
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Figure 6.- DSC block diagram. 
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Figure 11.- DSC and TGA purge gas and vacuum system. 









Figure l4.- Nitroglycerine infrared spectrum from Reference jk. 
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Figure 15.- Glyoxal infrared spectrum from Reference 35* 
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Figure 18.- Heat of reaction and burning rate versus pressure. 
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Figure 20.- Temperature jumps near surface (from ref. 10). 
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Figure 22. - TCA showing "knee" at same temperature as apparent 

endotherm DSC curve. 
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Figure 23.- Slow heating rate TGA curve which left carbonaceous 

residue. ^ 
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Figure 2k-. - Order of reaction, curve by method of reference 2/j. 
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Figure 25.- Computer printout near "knee" of curve 
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Figure 26- - Computer printout along second-order part of curve 
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SAMPLE CALCULATIONS 


1. Calculation of Activation Energy and Frequency Factor 
from DSC Data 

First, we will derive the expression to "be used. Since the rate 
constant is a function of temperature, it is different for each point 
on. our curve. However, the distance from the baseline to the heat of 
reaction curve (fig. 17) is proportional to the rate of heat evolution. 
It is therefore proportional to the rate constant, provided the 
Arrhenius plot is a straight line. If the values of the rate constant 
at Tj and Tg are kj_ and kg, we may write from equation (13) 
of the Arrhenius Theory section that 


E 


In k, = — — + constant 
1 RT- 


and 


E 


In kg = - + constant 

C Kip 

Subtracting (A-2) from (A-2) using Zn(x) - Zn(y) *= 


(A-1 ) ! 


(A-2) 


Zn ”, we obtain 

y 


Zn 


^ H ®1 T 2 _ 


(A-3) 


and since the distances d^ and dg at T^_ and Tg are proportional] 
to the rate constants 


Zn 


fa eQ _i“| 

dg - - % " Tgj 


(A-k) 



and solving for E 


-E = R 



U-5| 


Beginning at the point where the heat of reaction curve deviates from 
the Baseline;, make a plot of log d versus l/T. Use the extreme 
points of the straight line portion of this curve to substitute in 
(A-5). An example plot is given in figure 29 for the 50 mm Hg curve 
of figure 17« Then, using the values Z11 = 0.694, l/T^ = 0.002154; 

and Zn d 2 = 5-18, l/Tg = 0.002050 


-E = 1.986 « 47,000 ^ (A-6) 

(0.002154 - 0.002050) m0 ‘ L 


The straight line portion of the curve usually extends from the first 
recognizable deviation to the point near the top of the bell- shaped 
curve where the inflection of the curve begins to hend over the top 
of the peak. This procedure was repeated for several curves at each 
pressure to produce the data of table I. 

Notice that the calculation of the activation energy is indepen- 
dent of mass. If we know the mass of propellant which produced a 
given curve, we can calculate a frequency factor for the Arrhenius 
expression as follows: The distance from the baseline to the heat of 

reaction, curve of figure 17 is the heating rate in millicalaries per 
second produced by the mass of propellant present at that instant. 
Select a temperature along the bell- shaped portion of the curve which 
falls on the straight line portion of the Arrhenius plot as in 



figure 29’ At this temperature the heating rate in millicalories 
per second is given by 


q = m^f exp(- 


(A-T) 


where is the mass of propellant at the selected temperature I. i 
The fraction of the sample decomposed at temperature T is ^p/^total 


and we can write 


Am 

mj = m n 

•^total 


(A-8) 


where 


mQ = initial amount of mass 

A^t = area under the heat of reaction curve up to 
temperature T 

Apotal = area under heat of reaction curve 

The frequency factor is then given by 


Am _ 

1,0 Hot i E; 


*r 


(A -9) 


Using as an example the atmospheric pressure curve of figure 17 at 
T = 210° C (k&3° K), we have 


Sq = 1 mg 
_ 1 

^total 5 


H = UO ^ or 
2* gm 2rg 





Since 


k - A exp 


('til) 


A exp 


ax 

( E 1 _dt_ 

V " ■ x n 


Differentiating the logarithm of equation (A-13), we obtain 


Edt 

RT 2 


= d In 


(-sO- 


dm x 


(A-12) 


(A-13) 


(A-14- ) 


Integrating equation (A-l4 ) gives 


E 1 


/ dx\ 

('3t)- n 


A in X 


(A-13) 


Dividing equation (A- 10) by d in X and equation (A-ll) by A In Xj 
we get 

■(-*) 


d In 


Rdt 

RT £ ’d in X 


d In X 


- n 


(A-l 6) 


and 


(A-17) 


-HI) *»(-f) 

A In X A in X 
A plot of A(i/t)/a In X versus A ln(-dX/dt)j A In X should produce 
a straight line of slope -e/r and intercept n. 

In order to put this into a more usable form, let M equal the 
total number of moles in the reaction mixture and m a equal the moles 
of reactant A at time t. Than X = m^/M may be substituted in 
equation (A-13) and the same procedure followed to obtain 


In k = In Jyp 1 + - n Zn(m a ) 


(A-l8 ) 


±M 


f 

A H- if 


A In m a A In m a 


The weight or volume of the reactant may be measured to replace M 


and using 


where 


2% = .. dw 
dt w c dt 


W r = W c ~ w 


m 0 = initial, mole fraction of A 

w c = weight loss after reaction is completed 

w = total weight loss in time t 


LJ.YY 

*rr = reaction rate in time t 
dt 

Using equations (A-19)> (A-20), and (A-21), we obtain 


E a/ X \ a 

'b1t) , a Zn yitj 

A Zn w r A Zn w r 


A straight line should be obtained when A Zn(dw/dt )j A Zn w r is 
plotted versus A(l/T)j A in w r with intercept equal to the order of 
reaction n and slope equal to E/R. 



3- Order of Reaction Computer Program 




Performing the integration of the left side 



where 

P(X) = ~ - [-Ei(X)] 

and (A-29 ) 

ET 

can he 

(A-30) 

■tdiere 

\ = ”. 35^733 -, = 5 330657 

= 0.250621 h 2 = I.68153I1. 


The term Q-EiC-xJJ is the exponential integral. -E^(-X) 
approximated by^ 2 ^' 


- Ei (x) = X- 1 e -x/^ + *1* + ag\ 
\X2 + biX + bp / 



Therefore, 


~ln(- 


m - Hi- 


-*- ) for n = 


(n _ - (x - 0 -] 


for n ^ 1 


P(X) = K 
RT 


The computer program uses equation (A- 31 ) in an iteration scheme 
to find the parameters in the following sequence. The masses and 
temperatures at points from the experimental curve and a range of 
trial orders are the inputs to the program. 

Step 1. Select the first trial value of n. 

Step 2. Assume a trial value for E of 10,000. 

Step 3 . At the beginning temperature T^ of the selected intervs 


and at T 2 


% = or A x — 

1 TR 1 1 EPCX)! 


li ^ ^ ?(X) 2 or Ag = 


BP(X), 


For this value of n and tne starting value of E, A^_ and Aq axe 
computed. 

.Step Ij-. The program then iterates to find a value of E which 
will matte = Aq} that is, since • Aq = Ag 

™^1 _ ™2 ( A .-34 1 

EP(X) 1 EF(X) 2 



from which 


«g _ p(*)i 

K2 P(X) 2 


(A-35) 


The program balances this equation, using different E's in 



Step 5* These values of A and E are recorded for the assumed 
value of n which is held constant for each cycle. 

Step 6. Then, using the input mass values, compute (m/mo) for 
the entire curve. 

Step 7* Using these values, compute the errors 




(A-37) 
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and for n 1 


(*.) - -fN - _%) I ~ n + (n - DAI? /. if : 

\®o/ IBo \ Mo / * Tf \ RT/ 


(A-40) 


Step 8, At this time, we have a set of errors for each trial 
value of n of the form 



A linear fit is applied to the curve in the vicinity of where ej_ 
changes sign and this value of n where e-j. = 0 found. This is 
^inin' 5 M s and. the recomputed values of A and E are printed in 
the last row of each series of trial values of n as the calculated 
best fit. 

A discussion of the errors involved in using the first term, in 
the series approximation of the exponential integral is given by 
Nelson^ 2 ^. A listing of the Fortran statements for this program is 
given in appendix B. 

4- Computer Program for DSC Simulation 

The bell-shaped curves of figure 17 result from two effects. First 
heat is produced as a function of temperature at some exponential rate 
described by an Arrhenius expression. This can be described by a 


curve as follows: 




E = activation energy for heat release rate 
E = gas constant 
T - temperature 

All parameters in this expression are known from DSC measurements. The 
remaining portion of our desired expression must describe the mass 
loss rate, and it is obtained as follows: We have already measured 

the mass loss rate with the TGA system and evaluated it with the com- 
puter program described in the previous section. This program provides 
us with a means of expressing mass as a function of time of the proper 
order for our particular sample and heating rate range. Using the 
same assumed form of the rate equation (i.e., eq. (A-23)> follow the 
steps described in going to equation (A-28). Ry using an asymptotic 
expansion for -E^( -X), equation (A-29) can be written as 



If the first term of the series is defined as 


q(X) = X" 2 e" X (A-42) 


then K may be expressed as 


K = 


r “ 4 (I) 

RT 


(a-43) 


where 


r - 

q(X) 


(A-44 ) 



The term (l - r) is the relative error associated with using only the 
first term of the series. Nelson (^7) provides a plot of (1 - r) as a 
function of X which indicates that for X greater than 20. r is 
a slowly varying function approaching the value 1.0. Combining 
equations (A-42) and (A-43) gives 


K - rT 2 ~~ expf" 

ET V RT/ 


Evaluating equation (A-27) for our particular case with an apparent 
second-order process and taking m r = 0, we get 


/ m_ 

K = -2 - 1 
V m 


Setting this equal to equation (A-4-5) with r = 1 and solving for nij 
we obtain 


AR ’f E \ 
— r sxp ( - 1 I 
ET \ RT/ 


Combining this with our heat release expression and integrating with 
respect to time gives 


g< t > -/ 


o 1 + T 


EdV 

E*o \ 


H r f eK p(~ ^) dT 


(A-48) 
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where 


r = time in minutes 

Q(t) = the heat of the reaction at time t 
m Q = initial mass of sample 
T Q = temperature scan rate, °c/min 
Ed = activation energy for mass loss rate 
This equation, is solved on the digital computer by rewriting it as 


~ ■ 111 ~ 1 ■>- 1 " - ‘ - 1 " — : 

L J Mo + 1 „tJ 


V exp 


RjTi+Vf) 


dT - q(t) 


(A-49) 


where T has been replaced by + T q t . The integral is evaluated 
for different values of Ed between 20,000 and 60,000 cal/mole , using 
the TRAP subroutine shown as the last page in appendix B. The same 
iteration subroutine as used in the order of reaction program (i.e. , 
ITR-2) is then used to find which value of Ed in this range makes 
the equation equal aero. Ed is first found for a heating rate of 
T q = 20° C per minute. Using this value of Ed, the integral curve is 
plotted. This is considered as the reference curve. This same value 
of Ed is used again for T 0 = 10° and h0° C per minute. These 
curves are .plotted and compared to the reference curve. The resulting 
plots are shown in figure 31. 

The values used for the sample case are based on 1.0 mg of pro- 
pellant at atmospheric pressure. These values are: 

M 0 a 0.001 gm Tj = kk6° K T 0 = 20, 10, h0° K/min A = 10 2 3 

R = I.986 f = 10 21 E = 51,500 cal/mole H,. = k39 eal/gm 




PROGRAM I Tfr P t I NPL T ,QL TPLT ,1 APE5 = I NPLT ) 

INTEGER GRANGE 

REAL Ml ,M2 ,NC,.MR ,K1 ,X2 ,N « NF l N, NINC, NO, M, MMQE*M'MOC, NM IN 
C [MENS ION TC 30>,M( 30) ,MMOC< 30) ,MM3E(3G) ,N( 25) ,5UM(25) ,G(25 ,L) , 

1 M 25) , RES 10(25,1 ), SUB(1 ), ARR( 3 , 3) » B( 3 , 1 )» C C 25, 3) 

COMMON /RLK1/SLM,NN, CHANGE ,ERR,I 

CCNMCN/GLK?./ A0,A1, EC, El , K 1 ,K2, T DOT , E, A ,R , T 1 , T2, T , 1C0DE ,M,K 
NAMELI ST/NAM 1 /MO, MR, R ,TDDT » XNO , XNF I N ,X N) NC,XNN,K,M,T 
1 01 RFAC(5,1| ■ 

L FORMAT ( 5 1H I 

PRINT 1 
REAOt 5 ,NAM1 ) 

PRINT 3»M>J* MP, TOOT »R 

3 FCPMAT(1H0,3HMC=F7.3,5X,3HMR=F7 .3, 5X, 5HTD0T = F7.3, 5X, 2HR=F6.3//1 
PRINT E 

fl FORMAT I5X , 1 HN, 7X , 2F M 1 , 6X, 2 HM2, 10X.1HE, 17X , 1HA// ) 

MMIN=C. 

1=1 

NC1= XNO 
NFIN=XNF IN 
NI f'C=XM NC 
MN= XNN 
NCI) =NC 
J=1 
INC=l 

3 C C Tl=T( J) 

KI<=J + )NC 
T2 = Tt X Kl 
M1=M(J ) 

N.2 = M(KK) 

31 C CALL CAl.C(N( 1 ) , Mi , M2, MR, MO) _ 

I F ( I CODE. NF .0 ) GC TC 400 
PRINT S,N( I ),Ml,M2,E,A,SUM(I) ,ERR 
3 FORMAT (IX ,3F8 .3,4EI7 .8) 

GO TO 300 

400 FPINV 7, N ( I ), MI. ,MI2» ICODE 

7 F0RMATUX?3FR. 2.1RHNC SCLUT IGN, IC00E= I 2 > 

50 o IFIMD.GE.NFINJGO TO 600 
1=1 + 1 

M C I >=N ( I- 1 ) +N INC 
GC TC 300 

£CC CALL CHECK! ICH) 
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IF ( CHANGE, E G. 0 ) GC TC 200 

1= IO 

N0=NIIrll 

MNONINC/2. 

PvFI^=iNCI-»-L > 

NNt=<5 

1=1 

N< I )=N0 

330 CALL CALCtMU MP.MOJ 

I F C I CODE, NE . 0) GO TO 't00 
IF ( I .EC.NN1G0 TC 320 
1=1 + 1 

NIIHIH MNINC 
GO TC 230 
320 DO 30 n=l,NN 

win 1 = 1 . 

30 Gin,n=suMtm 

CALL LSCP0L(N»G,W.RESID,NN,SUB,1, ARK, B,2, C,25,3 » 
hM IN= l-Bll, 1H/E12, I) 

CALL CALCtNMlWMl ,V2 .('P.MGl 
. PRINT 5»NMIN,M1,M2,E»A»SUM(I) * ERR 
200 1=1 

PRINT C 
O FCRN AT tlHO) 

N{ I > = XNO 
N 1NC=XN INC 
NFIN=XNFI N 
NN=XNN 

IFl KK.EQ.KJGO 1 C 301 
J = J + 1 
GO TC 300 

3C1 IFl INC .FO . 3 )G0 TO 1C1 
I NC=I NC+1 
J= 1 

GC TC 300 
ENO 
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SUBROUTINE CALCl N ,M1, M2,PR » MO) 

Rfflt. M ,M2 , PC, «F,Kl,K?, M, NF IN» N INCt NO, N , ’IMO E ,.NM3C »NM IN 
COMPDW/RLKJ /SIM, NN, CHANGE, ERR, I 

COMMON'/ ELK?/ AC, A1 , Bi'lt 81, K 1 , K 2, T DOT , E » A ,R . T 1> TE.T.ICOOE ,M,K 
DIMENSION FMPCOOJ ,FMOE (30) , TI80), HI33 ) .SUM 125 ) 

EXTERNAL FCFX 

0 NT A AP,A1 ,BC»B1/ .250621,2.33733,1.681534, 3,3306 57/ 

Al=?. 334732 

DAT ft XLCW.XtP, El ,F2 ,MAX1,EELTX/103.»2F+09» 2*.0!301, 100000, SCO./ 
IFCWSIN. — 1 . 1 « LT. 0, 01) GO TO 10 

K 1 = 11 . / 1 N -l.n*U(Hl-MK)/MOJ**SJ, )-( I.-MR/MOI A *C 1 ,-N I ) 

K ?= { 1 . / ( M -1. U*< (( N2-MRJ/H0I-W-I1.-N )-(l.-MP./MI?»*#tl.-N II 

pc tt ?o 

1C K 1=-AL00( IK1-NRJ/1 N0-MR.I > 

K2 =— ALCC ( ( >/ INC-MR) ) 

2 n - CALL (E » XL 0> » XUP .DELTX, FGFX.El ,52 , MAXI , I CODE 1 

IT< ICf'CE.NE .O » GO TO 400 
>1=IE /( R^Tll 
X 2 - 0/ ( F *T 2 > 

1 c ( 48 SI X2 I . G T . 741 . 1 PRINT 1 ,E,R,T2,X2.N 
IF IABStX?) .GT.741. ) STOP 

1 FT u AI ( / /5 FI 6.5 //) 

E l 2— EXP { - X2 }/X2*< AD+M*X2+X2*X2J /IBO+B1* X2 + X2*X2) 

?= LK2* F+-TDCT )/ I E* i EXF (-X2 )/X2-E I2i ) 

IF(A6S(N -1-i.GT.O, CUGO TO 15 

F»C1-IPC-MR}/MC 

on JC^O L=1,K 

F AC.2 - FXP [ ( - A*T ( L MT { L >*£x P t~E/ f R*T ( L J M J ✓< T0UT*£ /R > ) 

MN|0C( U=MI</N0nFACI*F AC2 
1QH0 MMOF<L !=•-!(, L )/*«□ 

GC 1 7G ?5 

15 LACl-f (M0-MF J/'1CiJ**(l.-N ) 

OP 30PC 1=1 ,K 

FAC2=( A*T( L (*T!L J *EXP(-E/< R*T< DIM /( ( F/R> * TOOT) 

TEST-F ACl-M ft -1.MF6C2 
I F( TEST . GE . C . I GO TO 9C0 
4* MQC1L > =FR/ >C 
GO TO ?O0^ 

500 MMOCLL J-'TO/MO+l FACHIN -Ul*FHZ2S**i l./C 1.-N ») 

3CCC MMOFT L) ~'W ;L ) /V C 
?5 S 0 M ( J ;)= 0 . 

ERR-C. 

DO 2 COO L=1*K 
AOC=N , Mr FIl J-W0C1 L> 

Sq=Aon**2 
SIJM ( I) = SUM I I H- A HO 
2 CO C EP F=FRF+SO 

FRR=SOkT(FRR/<) 

4nn pETUFN 
EMO 
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FUNCTION FOFX(X) 

PFAL «2 ,K1 

CC M HON/rILK27AC,Al,ilO, 31 ,Kl , K2 ,TD(JT* ErA.K.Tl ,12, T , I CODE 
D I MENS ION T (30 I 
Xl=X/< K*T1) 

X 2= X / ( R *T 2 } 

FI l=FXP(-xn/Xl*(AC4Al*Xl+Xl*Xl )/ ( 8C+iil*X 1+X1*X1) 

F 1 2= EXP (-X? )/X2*< AO+A 1*X2+ X2+X2I /( 30+3 1* X2 + X2+X2) 

? = ( K2* F*TOCT)/(Xv {FXP (-X2 I/X2-EI2I ) 

FnFX=(Z*X/(R«TOOT) *(F XP(-Xl) /Xl-EI 1))-K1 
70C RETLPK . 

END 


5USPCUTINF CHECK t X C H > 

CHECKS FCR CP AKGE CF SIGN IN SOI 

INTEGER CFANGF 
0I1FKSTCK St'W(25I 
CCf'VCN/RLKl/SiJM, MM. CHANGE* ERR, I 
CH ANGF=C 

IFISUfM 1 ) .GE.O. )G0 TO 1 
GO TC 2 

1 DO 10 1=2 »NN 

IF <SU M m cGE.O.IGC TO 10 
CHANGF= 1 
TCH= I 

go rn ioo 
10 CCNTINUF 
GO TO ICO 

2 CO 20 1=2, NN 

I F ( SliM 1 1 J • i.1 » 0 « J GC TD 20 
CFAMGt= 1 
I CM=t 

GO TO IPn 
20 CCKTIMJE 
IOC RE Tliw K ‘ 

END 
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SUBROUTINE ITR2 


LANGUAGE: FORTRAN 

PURPOSE: Given F{X) = 0, to find a value for X within 

a given epsilon of relative error in a given 
interval (a,b). 

USE: CALL ITR2 (X, A, B, DELTX, FOFX, El, E2 , MAXI, ICODE) 

X The root. 

A The lower hound on X. This value is used by 
ITR2 as an initial guess. 

B The upper bound on X. This value is used by 
ITR2 as a final guess if the entire interval 
is scanned. 

DELTX AX, the size of the scanning interval. 

FOFX The name of a function subprogram to evaluate 
F(X). 

El Relative error criterion. 

E2 Absolute error criterion. 

MAXI A maximum iteration count supplied by the user. 

ICODE An integer supplied by ITR2 as an error code. 

This code should be tested by the user on re- 
turn to the calling program. 

ICODE = 0, normal return 

ICODE =1, maximum iterations are exceeded 
ICODE = 2, DELTX = 0, or negative 
ICODE = 3, a root cannot be found within 
t,he given hounds 
ICODE = 1< , A > B 

RESTRICTIONS : Make A<:B, AX positive. A function subprogram 

with a single argument X must be written by the' 
user to evaluate F(X). The name of this subprogram, 
FOFX 8 must Appear in an EXTERNAL statement of the 
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ITR2 


METHOD: 


ACCURACY; 


calling program. 

The given function F(X) is evaluated at a 
given starting point, a, and at intervals of 
a specified AX thereafter, up to and including 
a specified end point, b. A change of sign of 
the function across a AX interval indicates a 
possible root in that interval. The interval 
is then halved successively toward F(X) = 0 
until the prescribed accuracy is satisfied. 

The given function F{X) is evaluated once 
for each halving step. 

If the given function is expected to, have 
more than one root between the prescribed 
starting and end points, it is suggested 
that a sufficiently small AX be given 
such that no more than one root be present 
within a AX interval. A normal return is 
given upon the location of the first root 
from the starting point,- a. Additional 
roots mu ct be located by new entries into 
the subroutine using a new starting point, 
p, which is Just beyond the previous root. 


The iteration process is continued until 
either of two convergence criteria are 
satisfied. These criteria are: 


1 . 


if f f> Ci 

*i - *1-1 
Xi 


< € 


1 


REFERENCE: 

STORAGE: 


2, if X t < 

|*i - x i-i| < e 2 

T. i. Scarborough, NUMERICAL MATHEMATICAL 
TSIS, Fourth Edition. 

260 g locations 

NASA, LRG, Terry A. Street sr 


SOURCE: 



IDENTlFICATIOHi TRAP 


PURPOSE: 


RESTRICTIONS: 


USAGE: 


METHOD: 


ORIGmrOR: 


TRAP calculates the running Integral of the curve 

y * f(x). 


Y and X can "be input in table form or Y must be 
calculated at ML discrete points before entering TRAP. 

The subroutine has been ccmpilecL with a variable 
DIMENSION statement. The following must be iimensioned 
In the calling program: 

X(XL), Y(XL), SUM(IL). 


GAIL TRAP (Y, X, ML, SUM, SUMl) where 

Y = the name of the dependent variable. 

X = the name of the independent variable. 

ML « the number of points on the curve or in. the 
table. 

SUM = the location of a vector of order IL which 
contains the running integral of the curve 
y = f(x). SUM(XL) is the total integral. 

SUMl = SUM(l): The initial value of the integrated 

curve. 


TRAP uses the trapezoidal rule of Integration for a 
ve able AS. 


Raye C. Mathis 



